Soliton Approach to Spin-Peierls Antiferromagnets: Large-Scale Numerical Results 



00 



(N 



I 



I 

o 



00 

m 
in 
o 

00 
On 



o 
o 



X 



Erik S0rensen°, Ian Affleck'', David Augier" and Didier Poilblanc" 
Lahoratoire de Physique Quantique & UMR CNRS 5626, Universite Paul Sabatier, 31062 Toulouse, Franee 
''Department of Physics and Astronomy, and Canadian Institute for Advanced Research, University of British Columbia, 

Vancouver, BC, V6T IZl, Canada 
(February 1, 2008) 

A simple intuitive picture of spin-Peierls antiferromagnets arises from regarding the elementary 
excitations as S=l/2 solitons. In a strictly one-dimensional system these excitations are assumed 
not to form bound-states and to be repelled by impurities. Couplings to the three-dimensional 
lattice are assumed to produce an effective confining potential which binds solitons to antisolitons 
and to impurities, with the number of bound-states increasing as the interchain coupling goes to 
0. We investigate these various assumptions numerically in a phononless model where spontaneous 
dimerization arises from frustration and the interchain coupling is treated in mean field theory. 



The discovery of a spontaneous lattice dimerization in- 
duced by quasi one dimensional antiferromagnetic inter- 
actions, the spin-Peierls effect, in CuGeOa has sparked 
renewed experimental [1-10] and theoretical interest in 
this field. Subjects of current interest include the nature 
of the excitations in the dimerized phase [11-15], the ef- 
fect of impurities [16-18] and the effect of a magnetic 
field [19-21]. A simple intuitive picture of these phenom- 
ena is provided by the soliton model, first introduced 
in the context of frustrated spin chains by Shastry and 
Sutherland [22] and discussed for example by Khomskii 
and collaborators [23] as well as others [13,24-27]. This 
model is based on the assumption that all interchain in- 
teractions, both magnetic and elastic, are weak. While 
this assumption may not be very good for CuGeOa, it 
is possible that other spin-Peierls systems may be found 
which are more one-dimensional. 

In the absence of any explicit dimerization a com- 
pletely decoupled single S* = 1/2 chain can still have a 
spontaneously dimerized ground-state due to spin-phonon 
interactions or due to frustrating next nearest neighbor 
interactions. The ground-state is two-fold degenerate, 
corresponding to the two possible dimerizations, A and 
B, and the fundamental excitations are expected to be 
massive S=l/2 solitons, s, and antisolitons, s, with a gap 
Asoi, living on different sublattices, separating the two 
different dimerizations. It is sometimes assumed that no 
soliton-antisoliton bound-states form in this model and 
that solitons are repelled by non-magnetic impurities, i.e. 
by the ends of an open chain, as we demonstrate below. 
Neither of these assumptions is obvious and indeed some 
approaches make assumptions in contradiction to these 
ones [16]. 

It follows from general principles that spontaneous 
dimerization disappears at any non-zero T, due to a fi- 
nite density of solitons, for a strictly one-dimensional 
(ID) system. Interchain interactions, either magnetic or 
elastic, can fundamentally change the situation. Spon- 
taneous dimerization is now robust against the appear- 
ance of solitons. A soliton-antisoliton pair on a particular 



chain separated by some distance r leaves a region of the 
chain in the wrong ground-state relative to the neighbor- 
ing chains. The consequence is an energy cost Ar where 
A is determined by the interchain interactions. Such a 
linear potential is confining; the soliton and antisoliton 
cannot escape from each other and behave analogously to 
quark anti-quark pairs. Associated with this stabilization 
of the spontaneously dimerized phase is a drastic change 
in the excitation spectrum. As proposed by one of us [25] , 
low-lying excitations corresponding to soliton-antisoliton 
bound-states now have to occur. These excitations will 
have spin S* = 1 or and energy E = 2Asoi+Eii where Et, 
is the binding energy. There is no soliton continuum due 
to the confining nature of the interaction (i.e. the fact 
that it grows without limit at large separation), however 
a. ss — ss continuum can occur beginning at the lowest ly- 
ing S = 2 states. At this energy Ei, exceeds 2Asoi and it 
becomes energetically favorable to create a new ss pair. 
Furthermore, a soliton is bound by a linear potential to 
each non-magnetic defect. It is convenient to model the 
linear potential using a mean field treatment of interchain 
couplings. In this approach the interchain couplings sim- 
ply provide a term in the single-chain Hamiltonian which 
favors one or the other of the two dimerized states. It is 
important to realize that the soliton picture is only ex- 
pected to be useful if the interchain couplings are rela- 
tively small so that the solitons, although confined, still 
behave in a quasi-independent fashion at short distances. 
If the interchain couplings are large then the soliton and 
antisoliton never get far apart compared to their intrin- 
sic size (Compton wavelength) and behave effectively as 
a single well-defined magnon. Consequently, the soliton 
picture ceases to have much utility. 

In this paper we make an extensive numerical investi- 
gation of this soliton picture using the simplest possible 
model: 

H = J^[(l + ^(-1)^)S, • S,+i + J2S, • S,+i], (1) 

i 

a spin only Hamiltonian with spontaneous dimerization 
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arising from a second nearest neighbor interaction, J2, 
and explicit dimerization (representing the effect of the 
neighboring chains) due to an alternating nearest neigh- 
bor interaction, S. For (5 = this model exhibits a 
second order phase-transition to a dimerized phase at 
J2c — 0.2411 [28]. As a representative of the dimerized 
phase we shall take the vicinity of the Majumdar-Ghosh 
(MG) model [29], J2 = 5; since at this point the corre- 
lation length is vanishing and we therefore expect only 
minimal finite-size corrections. We have performed ex- 
act diagonalizations of up to 32 sites as well as density 
matrix renormalization group (DMRG) calculations with 
TO = 128 states. Our main results are as follows: When 
i5 = we find that there are no low energy bound-states 
near zero crystal momentum for J2 < 5- For J2 > 
the solitons behave as free massive particles [30]; we ex- 
plicitly calculate the gap and dispersion of the solitons 
and show numerically that the soliton is repelled by the 
ends of an open chain, behaving like a massive particle 
in a box. When 6 ^ we demonstrate that a ladder of 
S" = 1,0 soliton anti-soliton bound-states is formed, in- 
creasing in number with decreasing 5 and giving rise to 
a range of well-defined peaks in the dynamical structure 
factor close to g = tt. In this case, an isolated soliton will 
bind to one of the chain-ends. We begin by discussing 
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figure 1: a) The single soliton dispersion relation at the 
MG point for a 23 site chain. The solid line indicates 
the Shastry-Sutherland variational estimate uJsoi{k) = 
(J/2)(5/4-hcos2A;). b) < S'f > as a function of i. DMRG 
results are shown for L=101, with S^^^ = 1/2, m = 128. 

our results for 5 = 0. 

Odd length systems. The behavior of an isolated soli- 
ton can be studied by considering odd-length systems for 
which the ground-state has Stot = 1/2. In Fig. la we 
show results for the lowest-lying S^^^ = 1/2,3/2 states 
(open and full circles) of a 23 site chain at the MG point. 
The ground-state energy per spin of an even length sys- 




figure 2: The soliton gap as a function of J2. 



tern, for J2 = 1/2 and open boundary conditions, is 
—3 J/8 and we use this value to define zero energy for 
the odd-site system. A well-defined S=l/2 mode corre- 
sponding to the dispersion of the soliton is clearly visible 
around k = 7r/2. Approximating the soliton as a single 
unpaired spin between the two dimcr ground-states gives 
a rigorous upper bound on the soliton dispersion rela- 
tion [22,31] oi E = (J/2)[5/4 cos2fc]. This is shown 
as the solid line in Fig. 1 and agrees very well with the 
numerical data. Including 3 and 5 spin structures in the 
variational soliton wave-function reduces [32] the rigorous 
upper bound on the soliton gap, A^/J (at k = 7r/2) from 
0.125 to 0.11701 in remarkably good agreement with our 
best DMRG estimate of 0.1170(2) [30]. In Fig. lb < Sf > 
is shown for a L = 101 site system in the Sf^i = 1/2 
ground-state. Clearly the soliton is repelled by the open 
ends and enters approximately a particle in a box state < 
Sf >~ const-|-(— 1)* s\T?{-Ki/{L+l) indicated by the solid 
circles. We have calculated the soliton gap defined as the 
A,„i =\unL^^E{L+l)-{E{L) + E{L + 2))/2 {L even) 
as a function of J2. Our DMRG results are shown in 
Fig. 2. As J2 is increased from J2c the soliton gap should 
increase exponentially [25]: A^oi = exp(— fe/(J2 — J2c))- 
The numerical data seems largely consistent with such a 
behavior as can be seen in the inset of Fig. 2 although 
Aso/ quickly becomes too small for a reliable determina- 
tion. 

Even length systems. DMRG calculations can be per- 
formed using spin-inversion as a symmetry in which case 
even and odd multiplcts can be distinguished. Perform- 
ing such calculations for J2 = 1/2 we find that the gap 
to the lowest-lying triplet and singlet excitations are de- 
generate in the thermodynamic limit Aqo/J = Aoi/,/ = 
0.2340(2), precisely twice Ago;. This degeneracy per- 
sists for J2 < 1/2 and we conclude that there are no 
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low energy bound-states near zero crystal momentum, 
although such states could possibly occur for J2 > 1/2. 
On the other hand, at the MG point, Gaspers and Mag- 
nus (CM) [33] have shown that there are exact singlet and 
triplet bound-state at q = tt/2, degenerate with energy 
E/J = 1. The triplet state saturates the total weight 
of the dynamical structure factor which is a single delta 
peak at {q = 'k/2,E = J). The lowest lying S = 2 
state at q = 7r/2 marking the on-set of the continuum 
has E ~ 1.2 J and is clearly separated from the bound- 
state. Exact diagonalization results are consistent with 
the occurrence of bound-states for a small range of mo- 
menta close to (7 = 7r/2. It is interesting to compare this 
with the predictions of the sine-Gordon field theory, ex- 
pected to be valid near J2c- This relativistic field theory, 
for coupling constant, « Stt has no boimd-states. By 
Lorentz invariance, if there are no bound-states at zero 
momentum there cannot be any at finite momentum ei- 
ther. However, a resonance with a finite life-time would 
become longer-lived at finite momentum due to relativis- 
tic time dilation. Glearly the non Lorentz invariance at 
the MG point is important in allowing for bound-states 
only at finite momentum. 
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figure 3: a) The on-site magnetization for an even length 
chain, L = 100, at the MG point, in the two soliton 
state S^gf = 1 and parity P = —1. DMRG results with 
m = 128 for diflferent values of the explicit dimerization 
S is shown, m = 128 was used in the calculation. The 
solid line, for ^ = 0, indicates the uniform part of the 
magnetization, f) The on-site magnetization for an odd 
length chain, L = 51, at the MG point, with S^^^ = 1/2. 
DMRG results with m = 128 for diflferent values of the 
explicit dimerization S is shown. The chain begins with 
a weak link. 

Explicitly dim,erized systems. We next investigate the 
model with an alternating interaction, 5^0, added. 
Prom a numerical perspective two effects complicate such 
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figure 4: a) The lowest lying triplet (solid circles), singlets 
(open circles) and quintuplets (crosses) for J2 = 0.5, 5 = 
0.05, L = 28, as a function of fc = 2q. Results for L = 32, 
k = are shown to the left. The inset shows the same 
spectrum from J2 = 0.45, S = 0.10, L = 20. b) The 
dynamical structure factor S^^{q,ui) for J2 = 0.5, 6 = 
0.05, L = 28 and a broadening of e = 0.04J. The solid, 
dashed and dotted lines indicate the 3 triplet branches. 

an investigation; The exponentially diverging correlation 
length ^ as J2C is approached and the size, l^s, of a ss 
bound-state, diverging as (5 ^ 0. Note that Iss increases 
with the energy of the bound-state. We need to have 
L 3> ^, Igs- In Fig. 3a we show DMRG results for the on- 
site magnetization < Sf > in the lowest-lying two soliton 
state J2 = 1/2, S^^t = l, P = -1, L = 100 for various dif- 
ferent values of 6. For J two well-defined peaks can 
be seen in the uniform part of < 5? >, separated hy L/2. 
As 6 is increased the formation of a soliton bound-state is 
clearly visible and the excitation becomes more and more 
magnon like. In all cases the chain starts and ends with a 
strong link. Ghains starting and ending with a weak link 
will have edge like excitations [17]. In Fig. 3f we show 
DMRG results for < Sf > at the MG point for different 
values of 8 as calculated in the S^^^ — 1/2 ground-state. 
In this case L = 51 is odd and the first link of the chain is 
weak. Glearly the soliton binds to the end of the chain, 
in agreement with other recent results [17]. However, as 
5 is decreased the maximum in < Sf > moves further 
and further inside the chain, only for the largest value of 
5 is the maximum at the first chain site, as has also been 
foimd in the impurity susceptibility [17]. From these re- 
sults we can roughly estimate the size of a bound-state, 
Iss, as a function of 5. We see that we can only hope to 
determine the excited states if J2 ~ 1/2, ii 6 > 0.05 due 
to finite-size effects. 

We now examine the excitation spectrum for 6^0. 
In Fig. 4a we show the lowest lying triplet, singlet and 
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quintuplet excitations vs. k = 2q for a system with 
L = 28, J2 = 1/2,(5 = 0.05. Three well defined triplet 
branches are clearly visible below the continuum (the 
quintuplet S = 2 states) as well as two singlet branches. 
The third and highest lying singlet appears to remain 
marginally below the continuum at fc = as can be seen 
for the L = 32 results shown to the left of the panel. 
The inset in Fig. 4a shows the same spectrum but for 
L = 20, J2 = 0.45,5 = 0.1. In this case finite size cor- 
rections arc significantly smaller and we see that only 
two triplet and two singlet branches are below the con- 
tinuum. We take these results as clear evidence that 
the number of bound-states increases as 5 — > 0. At 
still larger 6 = 0.2, J2 = 0.35 only two triplets and a 
one singlet is found [15]. The point J2 = 0.45,5 = 0.1 
is along the disorder line 6=1 — 2J2 where the ex- 
cited states of Gaspers and Magnus [33] remain exact, 
as noted by Brehmer et al [13], however, they are no 
longer degenerate. These two states are indicated by an 
arrow in the inset (CM). As was the case at the MG 
point the triplet state saturates the structure factor at 
k = TT along the disorder line, i.e. S{k,u!) is a single 
delta peak at the energy of the triplet. In Fig. 4b we 
show results for the dynamical structure factor S^^ {q, iv) 
for the L = 28 results in panel a. Note the change be- 
tween k and q from panel a to b, = mr/l, q = mr/14. 
Here, S^^iq,Lo) = E„ \{^n\SM\^o)\^S,iu - E„+Eo), 
where Eq (En) is (arc) the energy (ics) of the ground- 
state $0 (triplet states $„), S^iq) = Yjje^'''^' Sj/VZ is 
the Fourier transform of SI and 5^ is a Lorentzian of 
width e. The three triplet branches arc most clearly vis- 
ible around q = n. Neutron scattering experiments on 
CuGeOs have so far only seen evidence for a single triplet 
branch at 5 = tt. However, the weight of a peak due to 
an eventual second triplet branch should be much smaller 
and could have been missed. Secondly, it is possible that 
other compounds may yield more clear evidence for a lad- 
der of soliton bound-states which would consolidate the 
soliton picture. In conclusion we have demonstrated that 
well defined soliton excitations occur in frustrated spin 
chains. In the absence of interchain coupling (6 = 0) the 
solitons do not bind and are repelled by the chain ends. 
In the presence of interchain coupling (5 7^ 0) a number 
of stable bound-states occurs and isolated solitons are 
attracted to the chain ends. We thank IDRIS (Orsay) 
for allocation of CPU time on the G94 and C98 CRAY 
supercomputers. The research of lA is supported in part 
by NSERC of Canada. 
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